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Abstract 

The cross section for the diffractive deep-inelastic scattering process ep — > eXp is mea- 
sured, with the leading final state proton detected in the HI Forward Proton Spectrometer. 
The data sample covers the range xp < 0.1 in fractional proton longitudinal momentum 
loss, 0.1 < \t\ < 0.7 GeV 2 in squared four-momentum transfer at the proton vertex and 
4 < Q 2 < 700 GeV 2 in photon virtuality. The cross section is measured four-fold differ- 
entially in t, xp, Q 2 and /3 = x/xp, where x is the Bjorken scaling variable. The t and xp 
dependences are interpreted in terms of an effective pomeron trajectory and a sub-leading 
exchange. The data are compared to perturbative QCD predictions at next-to-leading order 
based on diffractive partem distribution functions previously extracted from complementary 
measurements of inclusive diffractive deep-inelastic scattering. The ratio of the diffractive 
to the inclusive ep cross section is studied as a function of Q 2 , f3 and xp. 
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1 Introduction 



Diffractive processes such as ep — > eXp have been studied extensively in deep-inelastic elec- 
trons-proton scattering (DIS) at the HERA collider [fTUTTl. and provide an essential input for 
the understanding of quantum chromodynamics (QCD) at high parton densities. The photon 
virtuality Q 2 supplies a hard scale, which allows the application of perturbative QCD. Diffrac- 
tive DIS events can be viewed as resulting from processes in which the photon probes a net 
colour singlet combination of exchanged partons. A hard scattering QCD collinear factorisa- 
tion theorem [ 18 1 allows the definition of 'diffractive parton distribution functions' (DPDFs) for 
a given scattered proton four-momentum. The dependence of diffractive DIS on x, the Bjorken 
scaling variable , and Q 2 can thus be treated in a manner similar to inclusive DIS, e.g. through 
the application of the DGLAP parton evolution equations |fT9] - l22l . 

Within Regge phenomenology, diffractive cross sections are described by the exchange of 
a pomeron (P) trajectory, as illustrated in figure [Q In previous measurements at HERA |21[6l 
diffractive DIS cross sections are interpreted in a combined framework, which applies the QCD 
factorisation theorem to the x and Q 2 dependence and uses a Regge inspired approach to express 
the dependence on xp, the fraction of the incident proton longitudinal momentum carried by the 
colour singlet exchange. In this framework the data at low xp are well described and DPDFs 
and a pomeron trajectory intercept are extracted. In order to describe the data at larger xp, 
it is necessary to include a sub-leading exchange trajectory (M), with an intercept which is 
consistent with the approximately degenerate trajectories associated with the exchange of p, to, 
a 2 and f 2 mesons. 



Figure 1 : Schematic illustration of the diffractive DIS process ep — >■ eXp and the kinematic 
variables used for its description in a model in which the pomeron (]P) and a sub-leading (H) 
trajectory are exchanged. 

In many previous analyses diffractive DIS events are selected on the basis of the presence 
of a large rapidity gap (LRG) between the leading proton and the remainder of the hadronic 
final state X [[3l|5]|. The main advantage of the LRG method is high acceptance for diffractive 

'in this paper "electron" is used to denote both electron and positron. 
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processes. A complementary way to study diffractive processes is by direct measurement of the 
outgoing proton. This is achieved in HI using the Forward Proton Spectrometer (FPS) [|4l[T6l 
|23l . Although the FPS detector has low acceptance, the FPS method of studying diffraction has 
several advantages. In contrast to the LRG case, the squared four- momentum transfer t at the 
proton vertex can be reconstructed. The FPS method selects events in which the proton scatters 
elastically, whereas the LRG method does not distinguish between the case where the scattered 
proton remains intact or where it dissociates into a system of low mass M Y . The FPS method 
also allows measurements up to higher values of xp than is possible with the LRG method, 
extending into regions where the sub-leading trajectory is the dominant exchange. The FPS and 
LRG methods provide means to investigate whether the hard scattering process characterised 
by the variables f3 = x/xp and Q 2 depends also on the variables xp, t and My associated with 
the proton vertex. According to the proton vertex factorisation hypothesis, the cross section can 
be written as the product of two factors, one characterising the hard interaction depending on (3 
and Q 2 , the other characterising the proton vertex depending on xp and t. 

In this paper, a measurement of the cross section for the diffractive DIS process ep — > eXp is 
presented, using HI FPS data with statistics increased by a factor 20 compared to the previous 
analysis 01. In addition the kinematic range of the FPS measurement is extended to higher 
Q 2 . The high statistics of the present data make the FPS results competitive in precision with 
the results of the LRG method. Reduced diffractive cross sections, a^ 4: \(3,Q 2 ,xp } t) and 
a^ 3 \(3,Q 2 J xp), are measured. These measurements are used to to extract the parameters 
of the pomeron trajectory and to quantify the sub-leading exchange contribution. The proton 
vertex factorisation hypothesis is tested. The cross section dependence on the hard scattering 
variables, (3 and Q 2 , is further studied. The ratio of diffractive to inclusive ep cross sections is 
measured as a function of Q 2 , f3 and xp. The data are compared with similar measurements of 
the ZEUS experiment HI 1 LI 1 211 . The data are also compared directly with the LRG measurement 
in order to test the compatibility between the two measurement techniques and to quantify 
the proton dissociation contribution in the LRG data. 



2 Experimental Technique 

The data used in this analysis correspond to an integrated luminosity of 156.6 pb -1 and were 
collected with the HI detector in e~p interactions (luminosity of 77.2 pb _1 ) and e + p interactions 
(luminosity of 79.4 pb _1 ) during the HERA II running period from 2005 to 2007. During this 
period the HERA collider was operated at electron and proton beam energies of E e = 27.6 GeV 
and E p = 920 GeV, respectively, corresponding to an ep centre of mass energy of y/s = 
319 GeV. 



2.1 HI detector 

A detailed description of the HI detector can be found elsewhere H24U25L Here, the components 
most relevant for the present measurement are described briefly. A right handed coordinate 
system is employed with the origin at the position of the nominal interaction point that has its 
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z-axis pointing in the proton beam, or forward, direction and x(y) pointing in the horizontal 
(vertical) direction. Transverse momenta are measured with respect to the beam axis. 

The Central Tracking Detector (CTD), with a polar angle coverage of 20° < 9 < 160°, 
is used to reconstruct the interaction vertex and to measure the momentum of charged parti- 
cles from the curvature of their trajectories in the 1.16 T field provided by a superconducting 
solenoid. 

Scattered electrons with polar angles in the range 154° < 9' e < 176° are measured in 
a lead / scintillating-fibre calorimeter, the SpaCal ll26ll . The energy resolution is a(E)/E « 
7% I ^AB[GeVjffil% and the energy scale uncertainty is 1%. A Backward Proportional Chamber 
(BPC) in front of the SpaCal is used to measure the electron polar angle with a precision of 
1 mrad. Hadrons are measured in the Spacal with an energy scale precision of 7%. 

The finely segmented Liquid Argon (LAr) sampling calorimeter H27U281I surrounds the track- 
ing system and covers the range in polar angle 4° < 9 < 154°. The LAr calorimeter is used 
to reconstruct the scattered electron in DIS processes at high Q 2 . The LAr calorimeter consists 
of an electromagnetic section with lead as absorber, and a hadronic section with steel as ab- 
sorber. Its total depth varies with 9 between 4.5 and 8 interaction lengths. Its energy resolution, 
determined in test beam measurements, is o(E)/E « 11%/ WE[GeV] © 1% for electrons and 
a(E)/E Ki 50%/^E[GeV] © 2% for hadrons. The absolute electromagnetic energy scale is 
known to 1% precision. 

The hadronic final state is reconstructed using an energy flow algorithm which combines 
charged particles measured in the CTD with information from the SpaCal and LAr calorime- 
ters [|29ll . The absolute hadronic energy scale is known with a precision of 4% for the measure- 
ments presented here. 

The luminosity is determined with a precision of 3.7% by detecting photons from the Bethe- 
Heitler process ep — > epj in a calorimeter located at z = — 103 m. 

The energy and scattering angle of the leading proton are obtained from track measure- 
ments in the FPS Ifi~6ll23l . Protons scattered through small angles are deflected by the proton 
beam-line magnets into a system of detectors placed within the proton beam pipe inside two 
movable stations, known as Roman Pots. Each Roman Pot station contains four planes of five 
scintillating fibres, which together measure two orthogonal coordinates in the (x, y) plane. The 
fibre coordinate planes are sandwiched between planes of scintillator tiles used for the trigger. 
The stations approach the beam horizontally from outside the proton ring and are positioned at 
z = 61 m and z = 80 m. The detectors are sensitive to scattered protons which lose less than 
10% of their energy in the ep interaction and which are scattered through angles below 1 mrad. 

The leading proton energy resolution is approximately 5 GeV, independent of energy within 
the measured range. The absolute energy scale uncertainty is 1 GeV. The effective resolution 
in the reconstruction of the transverse momentum components of the scattered proton with 
respect to the incident proton is determined to be ~ 50 MeV for p x and ~ 150 MeV for p y , 
dominated by the intrinsic transverse momentum spread of the proton beam at the interaction 
point. The corresponding t-resolution varies over the measured range from 0.06 GeV 2 at \t\ = 
0.1 GeV 2 to 0.17 GeV 2 at \t\ = 0.7 GeV 2 . The calibration of the FPS is performed using a 
sample of elastic ep — > ep°p photoproduction events with p° — > it + it~ decays by comparing the 
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variables reconstructed in the CTD with the values measured in the FPS. The scale uncertainties 
in the transverse momentum measurements are 10 MeV for p x and 30 MeV for p y . Further 
details of the analysis of the FPS resolution and scale uncertainties can be found in @). For a 
leading proton which passes through both FPS stations, the average overall track reconstruction 
efficiency is 48%. 

2.2 Event selection and kinematic reconstruction 

The events used in this analysis are triggered on the basis of a coincidence between the FPS trig- 
ger scintillator tile signals and an electromagnetic cluster signal in the SpaCal or LAr calorime- 
ter. The trigger efficiency is around 99% on average. 

Several selection criteria are applied to the data in order to select the DIS event sample 
and to suppress beam related backgrounds, photoproduction processes and events in which the 
incoming electron loses significant energy through QED radiation. The DIS selection criteria 
are summarised below. 

• The reconstructed z coordinate of the event vertex is required to lie within 35 cm(~ 
3cr ) of the mean position. At least one track originating from the interaction vertex and 
reconstructed in the CTD is required to have a transverse momentum above 0.1 GeV. 

• The the energy E' e and the polar angle Q' e of the scattered electron are determined from the 
SpaCal (LAr) cluster, linked to a reconstructed charged particle track in the BPC (CTD), 
and the interaction vertex reconstructed in the CTD. The electron candidate is required to 
satisfy either 154° < 0' e < 176° and E' e > 8 GeV in the Spacal calorimeter or 0' e < 154° 
and E' e > 10 GeV in the LAr calorimeter. 

• The quantity E — p z , calculated from the energies and longitudinal momenta of all recon- 
structed particles including the electron, is required to lie between 35 GeV and 70 GeV. 
For neutral current DIS events this quantity is expected to be twice the electron beam 
energy neglecting detector effects and QED radiation. 

The following requirements are applied to the leading proton measured in the FPS. 

• The measurement is restricted to the region where the FPS acceptance is high by requiring 
the horizontal and vertical projections of the transverse momentum to lie in the ranges 
—0.63 < p x < —0.27 GeV and \p y \ < 0.8 GeV, respectively, and the energy of the 
leading proton E' to be greater than 90% of the proton beam energy E p , where E p is the 
energy of the proton beam. 

• The quantity E + p z , summed over all reconstructed particles including the leading pro- 
ton, is required to be below 1900 GeV. For neutral current DIS events this quantity is 
expected to be twice the proton beam energy neglecting detector effects. This require- 
ment is applied to suppress cases where a DIS event reconstructed in the central detector 
coincides with background in the FPS, for example due to an off-momentum beam proton 
(beam halo). 
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The inclusive DIS variables Q 2 , x and the inelasticity y are reconstructed by combining in- 
formation from the scattered electron and the hadronic final state using the method introduced 
in 10: 

* = rf + *(i-*) ; Q^'-Bm • * = S- U) 

Here, y e and y d denote the values of y obtained from the scattered electron only ('electron 
method') and from the angles of the electron and the hadronic final state ('double angle method'), 
respectively 113011 . 

Variables specific to diffractive DIS are defined as 

%F ~ q-P ' /3 -2g-(P-P')' () 

with q, P and P' denoting the four-vectors of the exchanged virtual photon and the incoming 
and outgoing proton, respectively. The variable (3 can be interpreted as the fraction of the 
longitudinal momentum of the colourless exchange which is carried by the struck quark. The 
variable x P is reconstructed directly from the energy of the leading proton, such that 

x P = l- E' p /E p . (3) 

Two methods are used to reconstruct (5 in order to obtain the optimal resolution across the 
measured xp range. It is reconstructed as = x/xp in the range xp > 0.012. For xp < 0.012 
the hadronic final state is used for the reconstruction according to: 



Q 2 
+ 

The mass Mx of the hadronic system X is obtained from 



* = -c?TW x ' (4) 



M x = (E 2 -pl-p 2 y -p 2 z ) had -^, (5) 

Vh 

where the subscript 'had' represents a sum over all hadronic final state particles excluding the 
leading proton and y^ is the value of y reconstructed using only the hadronic final state OTTl . 
Including the factor y/yn leads to cancellations of several measurement inaccuracies. 

The squared four-momentum transfer t = (P — P') 2 is reconstructed using the transverse 
momentum p t of the leading proton measured with the FPS and the value of xjp as described 
above, such that 



t = t ■ ^— t ■ = - p p (6) 



x w >m % 

1~ 111111 
- Xp l-Xjp 

where |t m i n | is the minimum kinematically accessible value of \t\ and m p is the proton mass. 

In this measurement, the reconstructed \t \ is required to lie in the range 0.1 < \t\ < 0.7 GeV 2 
and x ip in the range x P < 0.1. The measurement is restricted to a 'medium' Q 2 region of 
4 < Q 2 < 110 GeV 2 , 0.03 < y < 0.7 and a 'high' Q 2 region of 120 < Q 2 < 700 GeV 2 , 0.03 < 
y < 0.8 for data with electron candidates reconstructed in the Spacal and LAr calorimeters, 
respectively. The final data sample contains about 68 200 events at medium Q 2 and about 400 
events at high Q 2 . 
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3 Monte Carlo Simulation and Corrections to the Data 



Monte Carlo simulations are used to correct the data for the effects of detector acceptances, 
inefficiencies, migrations between measurement intervals due to finite resolutions and QED 
radiation. The reaction ep — > eXp is simulated with the RAPGAP program [32] using the HI 
2006 DPDF Fit B set Q. The HI 2006 DPDF Fit A and HI 2006 DPDF Fit B parameterisations 
give a consistent description of diffractive inclusive DIS processes flU, but the HI 2006 DPDF 
Fit B predictions are in better agreement with the diffractive di-jet production cross sections 
measured in DIS [6J. Contributions from both leading (P) and sub-leading (M) exchanges 
are considered. Hadronisation is simulated using the Lund string model lT3~3l implemented 
within the PYTHIA program [34J. An additional p-meson contribution relevant for the low Mx 
domain is simulated using the DIFFVM generator 11351 . 

The background from photoproduction processes, where the electron is scattered into the 
backward beam pipe and a particle from the hadronic final state fakes the electron signature, is 
estimated using the PHOJET Monte Carlo model [|36l . This background is negligible except at 
the highest y values and is 3% at most. The proton dissociation background, where the leading 
proton originates from the decay of a higher mass state, is estimated using an implementation in 
RAPGAP of the proton dissociation model originally developed for the DIFFVM Monte Carlo 
generator. This background is negligible except at the highest x$> values, where it reaches 2%. 

The response of the HI detector is simulated in detail using the GEANT3 program [|37l and 
the events are passed through the same analysis chain as is used for the data. 

Background mainly arises from random coincidences of DIS events resulting in activity in 
the central detector with off-momentum beam protons originating from interactions of beam 
protons with the residual gas in the beam-pipe or with the beam collimators (beam-halo back- 
ground) giving a signal in the FPS. This contribution is estimated statistically by combining the 
quantity E + p z for all reconstructed particles in the central detector in DIS events (without 
the requirement of a track in the FPS) with the quantity E + p z for beam-halo protons from 
randomly triggered events. The resulting background distribution is normalised to the FPS DIS 
data distribution in the range E + p z > 1900 GeV where beam-halo background dominates. 
After the selection cut of E + p z < 1900 GeV the background amounts to 13% on average. The 
E + p z spectra for leading proton and beam-halo DIS events are shown in figure [2^. The back- 
ground is determined using the reconstructed E + p z distribution as a function of the variables 
x, Q 2 and t. A comparison of the FPS data after background subtraction and the RAPGAP 
simulation is presented in figure [2}) for the energy of the leading proton E' p and in figure [3] for 
the variables xp, p x , p y and The beam-halo background is subtracted from the data. The 
Monte Carlo simulation reproduces the data within the experimental systematic uncertainties 
(section HI). 

Cross sections are obtained at the Born level by applying corrections for QED radiative ef- 
fects to the measured values. These corrrections amount to about 10% and are obtained using 
the HERACLES [38] program within the RAPGAP event generator. The measured cross sec- 
tions are quoted at the bin centres in Q 2 , (3,xp and \t\. Bin centre corrections are applied for 
the influence of the finite bin sizes using pomeron and sub-leading exchange parameterisations 
in the framework of the HI 2006 DPDF Fit B flU for the Q 2 , (3 and xp dependences and the 
measured t dependences at each (Q 2 , (5, xp) value (section ISTTI) . 
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4 Systematic Uncertainties on the Measured Cross Sections 



Systematic uncertainties are considered from the following sources. 

• The uncertainties in the leading proton energy and in the horizontal and vertical projec- 
tions of the proton transverse momentum are 1 GeV, 10 MeV and 30 MeV, respectively 
(section I2TTT) . The corresponding average uncertainties on the a? and a? measure- 
ments are 2.5%, 4.8% and 1.8%. The dominant uncertainty originates from the FPS 
acceptance variation as a function of the leading proton transverse momentum in the hor- 
izontal projection. 

• The electromagnetic energy scale uncertainty implies an error of 1% on the E' e mea- 
surement, which leads to an average systematic error of 1.2% on the cr? measurements. 
Possible biases in the 9' e measurement in the SpaCal (LAr) calorimeter at the level of 
±1 mrad (±3 mrad) lead to an average systematic error of 2.5%. 

• The systematic uncertainties arising from the hadronic final state reconstruction are de- 
termined by varying the hadronic energy scale of the LAr calorimeter by ±4% and that 
of the Spacal calorimeter by ±7%. These sources lead to an uncertainty in the a® mea- 
surements of around 1%. 

• The model dependence of the acceptance and migration corrections is estimated by vary- 
ing the shapes of the distributions in the kinematic variables xjp, f3, Q 2 and t in the 
RAPGAP simulation within the constraints imposed on those distributions to describe 
the present data. The xjp distribution is reweighted by (l/xjp) ±0 05 , the (3 distribution by 
^±0.05 and ^ _ ^)=fo.o5 5 the q2 distribution by log(Q 2 ) ±0 2 . These sources result in an 
uncertainty in the a® measurements of 1%. Reweighting the t distribution by e ±l results 
in an uncertainty of 1.4% for the measured range of 0.1 < |£| < 0.7 GeV 2 . 

• The model dependence of the bin centre corrections for the reduced cross section is esti- 
mated by comparing the results obtained in the framework of the HI 2006 DPDF Fit B 
and HI 2006 DPDF Fit A parameterisations flU for the kinematic variables and Q 2 . 
The xp parameterisation is reweighted by (l/xjp) ±0 05 . The average uncertainty for the 
reduced cross section is around 1%. Reweighting the t distribution by e ±l results in bin 
centre correction uncertainties of 1.6% for the measurements. 

• The uncertainties related to the subtraction of background are at most 2% for proton 
dissociation, 3% for photoproduction and 3% for the proton beam-halo contribution (sec- 
tion |3]). 

• The systematic error related to the reconstruction of the event vertex is on average 1%, 
as evaluated by comparing the reconstruction efficiencies for the data and Monte Carlo 
simulation. 

• A normalisation uncertainty of 1% is attributed to the trigger efficiencies (section [2i2l) . 
evaluated using event samples obtained with independent triggers. 
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• The uncertainty in the FPS track reconstruction efficiency results in the normalisation 
uncertainty of 2%. 

• A further normalisation uncertainty of 3.7% arises from the luminosity measurement. 

• The extrapolation in t from the measured FPS range of 0.1 < \t\ < 0.7 GeV 2 to the region 
|tmin| < |*| < 1 GeV 2 covered by the LRG data [5 J results in an additional normalisation 
error of 4% for the data (section 1531) . 



The systematic errors shown in the figures are obtained by adding in quadrature all contributions 
except the normalisation uncertainty, leading to an average uncertainty of 8% for the data. The 
overall normalisation uncertainties are of 4.3% and 6% for the ay 3 ^ and ay 3 ^ measurements, 
respectively. 



5 Results 

5.1 The reduced cross section cr^^ 

The dependence of diffractive DIS on f3, Q 2 , xp and t is studied in terms of the reduced diffrac- 
tive cross section . This observable is related to the measured differential cross section 
by 



d(3dQ 2 dx F dt f3Q 4 



[l-y + ^j-a^i^Q^x^t). (7) 



The reduced cross section depends on the diffractive structure functions F®^ and F^^ ac- 
cording to 

To a good approximation the reduced cross section is equal to the diffractive structure function 
F^ 4 \f3,Q 2 ,xp,t) in the region of relatively low y values covered by the current analysis. 
Results for ai° (4) are obtained in three t ranges, 0.1 < |t| < 0.3 GeV 2 , 0.3 < \t\ < 0.5 GeV 2 
and 0.5 < \t\ < 0.7 GeV 2 , and are interpolated to the values \t\ = 0.2, 0.4, 0.6 GeV 2 using the 
measured t dependence at each xjp, (3 and Q 2 value. Only the high statistics medium Q 2 data 
are used to evaluate the four-dimensional distributions a^^\ 

The reduced cross section xjpa^^ is presented in tabled Figure |4] shows xp&r^ as a 
function of xp for different \t\, [3 and Q 2 values. At medium and large values, xpa?^ 
falls or is fiat as a function of xp. Qualitatively this behaviour is consistent with a dominant 
contribution of the pomeron exchange described in the Regge framework by a linear trajectory 
ocp(t) = ap(0) + a'jpt with an intercept ap(0) > 1 [39J. At low values xpa? rises with 
xp at the highest xp, which can be interpreted as a contribution from a sub-leading exchange 
(M) with an intercept ap(0) < 1. This observation is consistent with the previous HI FPS 
analysis @J. 
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5.2 Cross section dependence on xjp and t and extraction of the pomeron 
trajectory 



The structure function is obtained by correcting ov^ 4 ' ) for the small F^^ contribution 

using the prediction of HI 2006 DPDF Fit B given in 0. To describe the xp and t dependence 
quantitatively, the structure function F®^ is parameterised by the form 

F 2 DW = f P {x P , t)F P {/3, Q 2 ) + n M ■ fmixjp, t)F M ((3, Q 2 ) , (9) 



which assumes proton vertex factorisation of the xp and t dependences from those on f3 and Q 2 
for both the pomeron and any sub-leading exchange, with no interference between the two con- 
tributions. The xp and t dependences are parameterised using flux factors fp and fp motivated 
by Regge phenomenology, 

Mx^ t) = Ajp. (xj3)2ajp(tM ; t) = A M . ^^pi , do) 

assuming that the sub-leading exchanges have a linear trajectory, ap(t) = ap(0) + a' M t, as 
is also assumed for the pomeron. Following the convention of (3), the values of Ap and Am 
are chosen such that xp ■ f* mm fp p(xp,t) dt = 1 at xp = 0.003 with t cut = —1 GeV 2 . 

Fitting the form of equation [9] to the experimental F®^ data, the free parameters in the fit 
are the intercept and slope of the pomeron trajectory, ap(t) = ap(0) + a' P t, the exponential 
t-slope parameter Bp for xp — > 1, the pomeron structure function Fp(f3, Q 2 ) at each of the 
Q 2 ) values considered, and the single parameter rip describing the normalisation of the 
sub-leading exchange contribution. As in j3]-0, the structure function Fp((3, Q 2 ) for the sub- 
leading exchange in each (3 and Q 2 bin are taken from a parameterisation of the pion structure 
function [|40l . 

The behaviour of F®^ at large xp and low is sensitive to the parameters of the sub- 
leading exchange 0:^(0), a' M and Bp. They are taken to be the same as in the previous fits to 
the HI F 2 D data [fflEl in order to compare the normalisation parameters for the the sub-leading 
exchange contribution between the measurements. The intercept Or(0) = 0.50 of the sub- 
leading exchange is taken from 0. The parameters a' M = 0.3 GeV~ 2 and Bp = 1.6 GeV -2 
are obtained from a parameterisation of the previously published HI FPS data [|4]|. The model 
dependent uncertainty is determined by repeating the fit with the fixed parameters made free 
one after another in the fit. The fitted parameters of the sub-leading exchange are consistent 
with the values given above. The influence of neglecting the FjP ^ contribution to is also 
included in the model dependent uncertainty. The experimental systematic uncertainties on the 
free parameters are evaluated by repeating the fit after shifting the data points according to each 
individual uncertainty listed in section @] 

The fit to equation [9] provides a good description of the x P and t dependences of the data, 
with a minimum x 2 = 273 for 289 degrees of freedom, combining statistical and uncorrelated 
systematic errors. The data hence support the proton vertex factorisation hypothesis for both 
the pomeron and the sub-leading contribution as given by the fit. 

The results for the free parameters of the fit are summarised in table Q] The experimental 
uncertainty of the fit parameters is defined as the quadratic sum of the statistical and systematic 
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Parameter 


Value 


ajp(O) 
a'jp 
Bp 
n M 


1.10 ± 0.02 (exp.) ± 0.03 (model) 
0.04 ± 0.02 (exp.) + J™ (model) GeV~ 2 
5.73 ± 0.25 (exp.) + ° f Q (model) GeV~ 2 
[0.87 ± 0.10 (exp.)+ (model)] ■ 10" 3 



Table 1 : The central values of the Regge model parameters extracted from a fit to F 2 1 ' and their 
experimental and model uncertainties. The experimental uncertainty is defined as the quadratic 
sum of the statistical and systematic uncertainties. The model uncertainty is determined by 
varying the fixed parameters in the fit as explained in the text. 

uncertainties. The overall normalisation uncertainty of ov^ 4 ' ) contributes only to the experimen- 
tal uncertainty of the sub-leading exchange normalisation parameter rip. The result for ap(0) 
is compatible with that obtained from HI data previously measured using the LRG and FPS 
methods [4,5] and with the ZEUS measurements H12U131I . It is also consistent with the pomeron 
intercept describing soft hadronic scattering, ap(0) ~ 1.08 H391I4UI421 . 

In a Regge approach with a single linear exchanged pomeron trajectory, ap(t) = ap(Q) + 
a'jpt, the exponential t-slope parameter B of the diffractive cross section is expected to decrease 
logarithmically with increasing x p according to 

B = Bp — 2a'p\iaxp , (11) 

an effect which is often referred to as 'shrinkage' of the diffractive peak. The degree of shrink- 
age depends on the slope of the pomeron trajectory, a'p. The present FPS data favour a small 
value of a'p, as expected in perturbative models of the pomeron [|43ll44l . This result is incom- 
patible with the value of a'p ~ 0.25 GeV~ 2 obtained from soft hadron-hadron scattering at 
high energies [|39ll4Tll42l . Vector meson measurements at HERA have also resulted in smaller 
values of a'p, whether a hard scale is present ||45l447ll or not l|48l. The present FPS results for 
a'p and Bp are compatible with those obtained previously from the HI / ZEUS data using the 
FPS / LPS detectors I141 PT21 Although the value of Bp measured in the HI experiment is lower 
than that from the ZEUS data (7.1 ± 0.7 (stat.) ±^7 (syst.) GeV~ 2 ), the results are consistent 
within uncertainties. 

The result for the sub-leading exchange normalisation parameter is slightly smaller but 
agrees within experimental uncertainties with rip = [1.0 ± 0.2 (exp.)] • 10~ 3 extracted from the 
previously published HI FPS data [4J and rip = [1.4 ± 0.4 (exp.)] • 10~ 3 obtained from the HI 
2006 DPDF Fit B to the HI LRG data flU. The FPS and LRG measurements give consistent 
results as expected for a dominantly isosinglet sub-leading trajectory (tu or f 2 , rather than p or 
a 2 exchanges). The sub-leading exchange is important at low j3 and high xp, contributing up to 
50% of the cross section at the highest bin centre value of xp = 0.075. 

5.3 Test of proton vertex factorisation 

To test in more detail a possible breakdown of proton vertex factorisation the dependence of 
ap(0), a'p and Bp on Q 2 is studied by repeating the fit described above in three different ranges 
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of Q 2 . The results of the fits, shown in figure[5]and table[2l indicate no strong dependence on Q 2 . 
The experimental uncertainty is defined as the quadratic sum of the statistical and uncorrelated 
systematic uncertainties. In the fit procedure, the normalisation factor n F for the sub-leading 
exchange contribution is fixed to the central value, presented in table Q3 as it is found to be 
insensitive to Q 2 . 



Q 2 range of Fit (GeV 2 ) 


MO) 


a'jp (GeV~ 2 ) 


Bjp (GeV" 2 ) 


4 < Q 2 < 12 
12 < Q 2 < 36 
36 < Q 2 < 110 


1.088 ±0.012 (exp.) 
1.102 ±0.016 (exp.) 
1.139 ±0.022 (exp.) 


0.009 ± 0.031 (exp.) 
0.063 ± 0.041 (exp.) 
0.023 ± 0.026 (exp.) 


5.78 ±0.20 (exp.) 
5.75 ±0.30 (exp.) 
5.17 ±0.40 (exp.) 



Table 2: The central values of a' F and B F and their experimental uncertainties extracted from 
fits to performed in three different ranges of Q 2 . 



In order to quantify a possible breakdown of proton vertex factorization, the data are fitted 
using parameterisations of the form A + D ■ \n(Q 2 /l GeV 2 ) for a F (0) and Bp. The logarithmic 
derivatives of a F (0) and B F are found to be D(a F (0)) = 0.018 ± 0.013 (exp.) and D(B F ) = 
—0.20 ± 0.14 (exp.) GeV -2 , respectively. Given the experimental uncertainties, the values of 
the logarithmic derivatives are within 1.5er from zero and hence do not contradict to proton 
vertex factorisation. 

Assuming an exponential t-dependence of the cross section, da/dt oc e Bt , the slope para- 
meter B is measured as function of x F at fixed values of Q 2 and j3. The results are presented in 
figure [6] The results for B are compared with a parameterisation of the t-dependence from the 
fit to (section |5?2~1) as shown in figured The fit results are shown as curves of the form: 

B(x F , (3, Q 2 ) = [1 - w M (x F , (3, Q 2 )} [Bp - 2a' F In x F ) + w M (x F , (3, Q 2 ) [B R - 2a' R In x F ] , 

w F (x F , f3, Q 2 ) being the fraction of F^^ which is due to the sub-leading exchange A good 
description of the data over the full x F , Q 2 and j3 range is obtained. At low x F , the data are 
compatible with a constant slope parameter, B ~ 6 GeV -2 . No significant Q 2 or (3 dependence 
of the slope parameter B is observed for data points with x F < 0.025. The sub-leading ex- 
change contribution integrated over this kinematic range is 5%. A parameterisation of the data 
in this x F range with a constant slope parameter B gives x 2 — 89 for 75 data points, where 
the errors include the combined statistical and uncorrelated systematic uncertainties. Within 
uncertainties, the t dependence of the cross section in the pomeron dominated low x F region 
can therefore be factorised from the Q 2 and (3 dependences. 

Since no significant Q 2 or (3 dependence is observed, the slope parameter B is obtained 
by averaging over the Q 2 and (3 and compared with the result of a parameterisation of the 
t-dependence from the fit to F% . The result is shown as a function of x F in figure |7J The 
previously published HI FPS results [4] are also shown. A weak decrease of the B parameter 
value from 6 GeV -2 to less than 5 GeV -2 is observed towards larger values of x F >0.05, 
where the contribution from the sub-leading exchange is significant. This reduction of the 
slope parameter indicates that the size of the interaction region is reduced for M, exchange, as 
compared to IP exchange. 
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5.4 The reduced cross section cr^ 3 ) and comparison with other measure- 
ments 



The reduced cross section a^ 3 \/3, Q 2 ,xp), defined as the integral of <Jr^\/3, Q 2 ,xp,t) overt 
in the range |t m i n | < \t\ < 1 GeV 2 , is obtained by extrapolating the FPS data from the measured 
ranged < \t\ < 0.7 GeV 2 using the t-dependence at each (xp, 3. Q 2 ) value (section l5.3l) . This 
extrapolation factor, which amounts to a value of 1.8 with an uncertainty of 4%, depends only 
weakly on xp. 

In figure [8] the HI FPS measurements of Xpert? are presented. They are compared with 
those of the ZEUS collaboration, measured using their Leading Proton Spectrometer (LPS) lfi~2ll . 
The ZEUS data points are interpolated to the j3, Q 2 and xp values of this measurement using 
a parameterisation of the ZEUS DPDF SJ fit [15 J. The ratio of the HI FPS to ZEUS LPS data 
averaged over the measured kinematic range is 0.85 ± 0.01(stat.) ± 0.03(syst.)^Q ^(norm.), 
which is consistent with unity taking into account the normalisation uncertainties of 6% and 
for the HI FPS and ZEUS LPS data, respectively. Within the errors, there is no strong 
xp, (3 or Q 2 dependence of the ratio. The FPS data extend the kinematic range of the cross 
section measurement to higher Q 2 and low (3. 

The reduced cross section can be compared with HI measurements obtained using the 
LRG technique J51 after taking into account the slightly different cross section definitions in the 
two cases. The cross section ep — > eXY measured with the LRG data is defined to include pro- 
ton dissociation to any system Y with a mass in the range My < 1.6 GeV, whereas Y is defined 
to be a proton in the cross section measured with the FPS. The results on xp a? measured 
using the FPS and LRG methods are shown in figure|9]as a function Q 2 in bins of and xp, and 
in figure [10] as a function of in selected bins of Q 2 and xp. The kinematic range of the mea- 
surements is extended to higher xp. The experimental uncertainties of the two measurements 
are defined as the quadratic sum of the statistical and uncorrected systematic uncertainties. As 
can be seen in figure [TQl the present FPS measurement has a precision comparable to the mea- 
surement [5] obtained using the LRG method. The LRG results are interpolated to the Q 2 , (3 
and xp bin centre values of the FPS data using the parameterisation HI 2006 DPDF Fit B Q. 

Since the two data sets are statistically independent and the dominant sources of systematic 
errors are different, correlations between the uncertainties of the FPS and LRG data are neg- 
ligible. The ratio of the two measurements is formed for each (Q 2 , f3, xp) point in the range 
xp < 0.04, where LRG data are available. The dependence of this ratio on each kinematic 
variable is studied by taking statistically weighted averages over the other two variables. 

The ratio of the LRG to the FPS cross section is shown in figure [TT| as a function of Q 2 , 
(3 and xp. Within the uncorrected uncertainties of typically 6% per data point, there is no 
significant dependence on (3, Q 2 or xp. The ratio of overall normalisations, LRG / FPS, is 
1.18 ± 0.01(stat.) ± 0.06(uncor.syst.) ± 0.10(norm.), the dominant uncertainties arising from 
the normalisations of the FPS and LRG data. This result is in agreement within uncertainties 
with the value of 1.23 ± 0.03(stat.) ± 0.16(syst.) obtained from the previously published HI 
LRG and FPS data @]|. Combining the result of |4] with the present measurement leads to a 
more precise value of the cross section ratio: 

a(M Y < 1.6 GeV) 

-±-£- = 1.20 ± 0.11 exp. • (12) 

a[My = m p ) 
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where the experimental uncertainty is a combination of the statistical, uncorrelated systematic 
and normalisation uncertainties of the measurements. The result is consistent with the predic- 
tion of 1.15+q from the DIFFVM generator, where the total proton elastic and proton disso- 
ciation cross sections are taken to be equal for the central value and their ratio is varied in the 
range 1 : 2 to 2 : 1 for the uncertainties [1511351 . 

The good agreement, after accounting for proton dissociation, between the LRG and the 
FPS data confirms that the two measurement methods lead to compatible results, despite their 
very different systematics. The lack of any kinematic dependence of the ratio of the two cross 
sections shows, within the uncertainties, that proton dissociation with My < 1.6 GeV can be 
treated in a similar way to the elastic proton case. The result confirms that contributions from 
proton dissociation in the LRG measurement do not significantly alter the measured j3, Q 2 or 
xjp dependences and hence cannot have a large influence on the diffractive parton densities 
extracted from the LRG data up to the normalisation difference. 

5.5 Cross section dependence on Q 2 and (3 

The measured reduced cross sections xpa^ 3 ^ are presented in figures [T2lfT4l and in table [4] 
as a function of xp, (3 and Q 2 . The FPS data are compared with QCD predictions at next-to- 
leading order derived from the HI 2006 DPDF Fit B to the HI LRG cross sections JH, which 
include both the pomeron and a sub-leading exchange. The normalisation of the HI 2006 
DPDF Fit B predictions is reduced by a global factor 1.20 to correct for the contributions of 
proton dissociation processes to the LRG cross sections, as evaluated in section [541 

As can be seen in figure [121 the rise of the data at large xp is in agreement with a significant 
contribution from a sub-leading exchange. The reduced cross section a? shown in figure [T3l 
decreases with (3 over most of the kinematic range. However, it clearly rises as /3 — > 1 at low 
Q 2 and xp. Within the framework of DPDFs, this can be explained in terms of diffractive quark 
densities peaking at high fractional momentum at low Q 2 [|3l[5]|. 

The figure [ID shows the Q 2 dependence of oy^ 3 ' at fixed xp and (3. Positive scaling vi- 
olations (d /din Q 2 > 0) are observed throughout the kinematic range, except at the 
highest (3 = 0.56. This observation is consistent with previous HI measurements using the 
LRG method (3l[5l and implies a large gluonic component in the DPDFs. As can be seen from 
QCD predictions, the positive scaling violations may be attributed to the pomeron contribution 
even at the highest xp values, where the sub-leading exchange is largest. The Q 2 dependence 
is quantified by fitting the data at fixed xp and (3 to the form 

xpa?®(/3,Xp,Q 2 ) = a D (/3 1 xp) + b D (/3,xp)\n(Q 2 /l GeV 2 ) (13) 

such that b]j((3,Xp) is the derivative of the reduced cross section with respect to h\Q 2 multiplied 
by xp. This form is fitted to data points for which the xp bin centre values satisfy x P < 0.035 
and for which the (3 bin contains at least 3 data points. The sub-leading exchange contribution 
at xp = 0.035 is below 15%. The resulting logarithmic derivatives are shown in figure [TBI 
Although the logarithmic derivatives at different xp values cover different Q 2 regions, they are 
similar when viewed as a function of (3. This confirms the applicability of the proton vertex 
factorisation framework to the description of the current data. The FPS results are consistent 
with predictions derived from the HI 2006 DPDF Fit B to the HI LRG data also shown in figure 

ED 
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5.6 Comparison of the diffractive and inclusive DIS cross sections 

By analogy with hadronic scattering, the diffractive and the total cross sections in DIS can be 
related via the generalisation of the optical theorem to virtual photon scattering [|49l . Many 
models of low x DIS ll50T - [55l assume links between these quantities. Comparing the Q 2 and 
x dynamics of the diffractive with the inclusive cross section is therefore a powerful means of 
testing models and of comparing the properties of the DPDFs with their inclusive counterparts. 
A detailed comparison of the diffractive and inclusive cross section is performed in [|5]|. Fol- 
lowing 0, the evolution of the reduced diffractive cross section with Q 2 is compared with that 

.D(3) / a /~\2\ 

of the reduced inclusive DIS cross section ay by forming the quantity (1 — fi Y F °\x=p^Q?) a * 
fixed Q 2 , (3 and xp, using a parameterisation of the ay data from ll56ll . This quantity is equiva- 
lent to the ratio of diffractive to inclusive DIS cross section M\ daD(K >{ d ^ x / W ' Q2) /o]*^ X {W, Q 2 ) 
studied in lfTTl[T3l[T4l as a function of the YP centre of mass energy W and Q 2 in ranges of 
Mx- The ratio is shown in figure [16] as a function of (3 at fixed xp and Q 2 . 

The ratio of the diffractive to the inclusive cross section is approximately constant as a 
function of (3 at fixed Q 2 and xp except at high 0. As can be seen in figure [161 the decrease of 
the ratio towards high (3 is reproduced by QCD predictions based on diffractive and inclusive 
proton PDFs [[5l[56]|. The ratio also rises towards larger values of xp where the sub-leading 
exchange contribution to the diffractive cross section is not negligible. 

The ratio, shown in figure [FT] as a function of Q 2 at fixed xp and f3, depends only weakly 
on Q 2 for most (3 and xp values. In order to compare the Q 2 dependences of the diffractive 
and the inclusive cross sections quantitatively, the derivative b r of the ratio with respect to 
\nQ 2 is extracted through fits of the form a r ((3, xp) + b r (/3, xp) ■ ln(Q 2 /l GeV 2 ). To reduce 
the influence of the sub-leading exchange contribution, only data points with bin centres at 
xp < 0.035 are included in the analysis of b r . The resulting values of b r are shown in figure [T8l 
They are consistent with zero within 3a. At fixed (3, there is no significant dependence of the 
logarithmic derivative on xp. Whereas the reduced diffractive and inclusive cross sections are 
closely related to their respective quark densities, the logarithmic derivatives are approximately 
proportional to the relevant gluon densities in regions where the Q 2 evolution is dominated by 
the g — >■ qq splitting [57] . The compatibility of b r with zero thus implies that the ratio of the 
quark to the gluon density is similar in the diffractive and inclusive DIS when considered at 
the same low x value. As can be seen in figure [T8l QCD predictions based on proton PDFs 
extracted in diffractive and inclusive DIS (51[56]| reproduce a weak decrease of the logarithmic 
derivative towards larger (3. 



6 Summary 

A cross section measurement is presented for the diffractive deep-inelastic scattering process 
ep — > eXp. The results are obtained using high statistics data taken with the HI detector at 
HERA. In the process studied, the scattered proton carries at least 90% of the incoming proton 
momentum and is measured in the Forward Proton Spectrometer (FPS). The data cover the 
range xp < 0.1 in fractional proton longitudinal momentum loss, 0.1 < \t\ < 0.7 GeV 2 in 
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squared four-momentum transfer at the proton vertex, 4 < Q 2 < 700 GeV 2 in photon virtuality 
and 0.001 < (3 = x/xp < 1. The measurement is performed in the range of the inelasticity 
variable 0.03 < y < 0.7 for 4 < Q 2 < 110 GeV 2 and in the range 0.03 < y < 0.8 for 120 < 
Q 2 < 700 GeV 2 . The new HI FPS data are in good agreement with the previously published 
HI FPS results and are consistent within uncertainties with results of the ZEUS collaboration 
obtained with their Leading Proton Spectrometer. The new measurements extend the kinematic 
range to higher Q 2 values. 

The reduced diffractive cross section a?^\l3, Q 2 , x F , t) is measured. The xp and t de- 
pendences are described using a model which is motivated by Regge phenomenology, in which 
a leading pomeron and a sub-leading exchange contribute. The effective pomeron intercept 
describing the data is atp(0) = 1.10 ± 0.02 (exp.) ± 0.03 (model), which is compati- 
ble within uncertainties with the pomeron intercept measured in soft hadron-hadron scatter- 
ing. The slope of the pomeron trajectory a'p is consistent with zero and smaller than the 
value ~ 0.25 GeV -2 obtained from soft hadron-hadron scattering data. The t-dependence 
of the pomeron exchange is described by an exponential function with constant slope param- 
eter Bp = 5.73 ± 0.25 (exp.) *°'g (model) GeV" 2 . The measured values of the slope of 
the pomeron trajectory and the t-dependences are characteristic of diffractive hard scattering 
processes. The Q 2 dependence of the parameters ajp(0), a'p and Bp is studied. The logarith- 
mic derivatives are consistent with zero within 1.5er of the experimental uncertainties thereby 
supporting the proton vertex factorisation hypothesis. 

The data are also analysed in terms of the reduced diffractive cross section a^ 3 \ obtained 
by integrating over the range |t min | < \t\ < 1 GeV 2 . At fixed x F , a decrease of <t, D( - 3 ' ) 
with j3 is observed over most of the kinematic range, except at the lowest values of Q 2 and 
xp. The data display positive scaling violations except at the highest j3 value of 0.56. The size 
of the measured scaling violations implies a large gluonic component in the diffractive parton 
distributions, in agreement with previous observations. 

The FPS results are compared with those obtained from an earlier HI measurement using 
events selected on the basis of a large rapidity gap rather than a leading proton. This LRG 
measurement includes proton dissociation to states Y with masses M Y < 1.6 GeV. The FPS 
data extend the kinematic range to higher xp and thus constrain the sub-leading exchange 
contribution. The ratio of the LRG to the FPS cross section is 1.20 ± 0.11 (exp.). It is 
independent of Q 2 , (3 and xp within uncertainties, confirming that contributions from proton 
dissociation in the LRG measurement do not significantly alter the measured Q 2 , (5 or x P 
dependences. 

The ratio of the diffractive to the inclusive ep cross sections is measured as a function of 
Q 2 , (5 and xp. At fixed xp the ratio depends only weakly on Q 2 or on j3 except at the highest 
(3. QCD predictions based on diffractive and inclusive proton PDFs reproduce the behaviour of 
the ratio. This result implies that the ratio of quark to gluon distributions at low x is similar in 
the diffractive and inclusive processes. 
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Table 4: The reduced diffractive cross sections xpo r as a function of Q 2 ,/3 and xp values (columns 1 — 4). The statistical (S sta t), 
systematic (5 sys ) and total (5 to t) uncertainties are given in columns 6 to 8. The remaining columns give the changes of the cross sections 
due to a +la variation of the various systematic error sources described in section|4j the hadronic energy scale (5had)\ the electromagnetic 
energy scale (<5 e ; e ); the scattering angle of the electron (5e); the reweighting of the simulation in /3 (5p), xjp (S Xp ) and t (5 t ); the leading 
proton energy E p (Se p ), the proton transverse momentum components p x (S Px ) and p y (S Py ); the reweighting of the simulation in Q 2 (8q2); 
the background from beam halo, photoproduction and proton dissociation processes (<5b 9 „); the vertex reconstruction efficiency (5 vtx ) and 
the bin centre corrections (5& cc ). All uncertainties are given in per cent. The normalisation uncertainty of 6% is not included. The table 
continues on the next pages. 
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Figure 2: (a) The distribution of E + p z for FPS DIS events (histogram with error bars) and 
for random coincidences of DIS events reconstructed in the HI central detector with beam-halo 
protons giving a signal in the FPS (histogram with shaded bands). The systematic uncertainties 
on the beam-halo background are presented as shaded bands around the beam-halo histogram, 
(b) The distribution of the leading proton energy reconstructed in the FPS (histogram with error 
bars). The beam-halo background is subtracted from the data. The RAPGAP Monte Carlo 
simulation is shown as a histogram with shaded bands indicating the experimental systematic 
uncertainties. 
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Figure 3: The distributions of the variables (a) xjp, (b) p x , (c) p y and (d) \t\ reconstructed using 
the FPS (histogram with error bars). The beam-halo background is subtracted from the data. 
The RAPGAP Monte Carlo simulation is shown as a histogram. The experimental systematic 
uncertainties are presented as shaded bands around the Monte Carlo histogram. 
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Figure 4: The reduced diffractive cross section xjpar^(P, Q 2 , xjp, t) shown as a function of 
xjp for different values of t, (3 and Q 2 . The error bars indicate the statistical and systematic 
errors added in quadrature. The overall normalisation uncertainty of 4.3% is not shown. The 
solid curves represent the results of the phenomenological Regge fit to the data, including both 
the pomeron (iP) and a sub-leading (M) exchange. 
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Figure 5: Results for (a) ajp(O), (b) a'jp and (c) Bp obtained from a modified version of the 
Regge fit performed in three different ranges of Q 2 . The error bars correspond to the experi- 
mental uncertainties which are the statistical and uncorrelated systematic uncertainties added 
in quadrature. The white lines and shaded bands show the result and experimental uncertainty 
from the standard fit over the whole Q 2 range. 
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Figure 6: Results for the slope parameter B obtained from a fit of the form da/dt oc e Bt shown 
as a function of xjp for different values of (3 and Q 2 . The error bars indicate the statistical and 
systematic errors added in quadrature. The solid curves represent the results of the phenomeno- 
logical Regge fit to F 2 D< ~ 4 ' ) including both the pomeron (IP) and a sub-leading (M) exchange. 
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Figure 7: The slope parameter B obtained from a fit of the form da/dt oc e Bt shown as a 
function of x#>. The data are averaged over Q 2 and (3. The inner error bars represent the 
statistical errors. The outer error bars indicate the statistical and systematic errors added in 
quadrature. The solid curve represents the results of the phenomenological Regge fit to the 
data, including both the pomeron (iP) and a sub-leading (M) exchange. The dashed curve 
represents the prediction beyond the x p range used in the fit. The previously published HI FPS 
results JH are also shown (open circles). 
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Figure 8: The reduced diffractive cross section xjpo r ((3,Q 2 ,xip) for \t\ < 1 GeV 2 , shown 
as a function of xp for different values of j3 and Q 2 . The error bars indicate the statistical and 
systematic errors added in quadrature. The HI FPS data are compared with the ZEUS LPS 
results |[T2l interpolated to the FPS f3, Q 2 ,xjp values. The overall normalisation uncertainty of 
6% on the HI FPS data and the normalisation uncertainty of on the ZEUS LPS data are 
not shown. 
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Figure 9: The reduced diffractive cross section xpa r (/3,Q 2 ,£jp) for \t\ < 1 GeV 2 , shown 
as a function of Q 2 for different values of xp and 0. The error bars indicate the statistical and 
systematic errors added in quadrature. The HI FPS data are compared with the HI LRG results 
interpolated to the FPS f3, Q 2 , xp values [5]. The overall normalisation uncertainty of 6% on 
the FPS data and the normalisation uncertainty of 6.2% on the LRG data are not shown. 
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Figure 10: The reduced diffractive cross section xjp a r {/3, Q 2 , xjp) for \t\ < 1 GeV 2 , shown 
as a function of for selected values of xp and Q 2 . The inner error bars represent the statistical 
errors. The outer error bars indicate the statistical and systematic errors added in quadrature. 
The HI FPS data are compared with the HI LRG results [5| interpolated to the FPS (3,Q 2 ,X]p 
values. The solid curves represent HI 2006 DPDF Fit B to the LRG data. The overall normali- 
sation uncertainty of 6% on the FPS data and the normalisation uncertainty of 6.2% on the LRG 
data are not shown. 
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Figure 11: The ratio of the reduced diffractive cross section a r for My < 1.6 GeV and 
\t\ < 1 GeV 2 obtained using the HI LRG data to that for M Y = m p and \t\ < 1 GeV 2 , 
obtained from the present and previously published FPS data [4J. The results are shown as a 
function of (a) Q 2 , (b) (3 and (c) xp, after averaging over the other variables. The inner error 
bars represent the statistical errors. The outer error bars indicate the statistical and uncorrected 
systematic errors added in quadrature. The combined normalisation uncertainty of 8.5% is 
shown as a band in figure (c). The dashed line in figure (c) represents the result of a fit to the 
data in the region shown assuming no dependence on xjp. 
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Figure 12: The reduced diffractive cross section xpo r (/3, Q 2 , xjp) for \t\ < 1 GeV 2 , shown 
as a function of xp for different values of (3 and Q 2 - The error bars indicate the statistical 
and systematic errors added in quadrature. The overall normalisation uncertainty of 6% is not 
shown. The solid curves represent the results of the HI 2006 DPDF Fit B to the LRG data J51 
reduced by a global factor 1.20 to correct for the contributions of proton dissociation processes. 
The dashed curves indicate the contribution of pomeron exchange in this model. 
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Figure 13: The reduced diffractive cross section xpo r (/3, Q 2 , xjp) for \t\ < 1 GeV 2 , shown 
as a function of for different values of xjp and Q 2 . The error bars indicate the statistical 
and systematic errors added in quadrature. The overall normalisation uncertainty of 6% is not 
shown. The solid curves represent the results of the HI 2006 DPDF Fit B to the LRG data J51 
reduced by a global factor 1.20 to correct for the contributions of proton dissociation processes. 
The dashed curves indicate the contribution of pomeron exchange in this model. 
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Figure 14: The reduced diffractive cross section xjp o r (/?, Q 2 , £p) for \t\ < 1 GeV 2 , shown 
as a function of Q 2 for different values of xjp and /3. The error bars indicate the statistical 
and systematic errors added in quadrature. The overall normalisation uncertainty of 6% is not 
shown. The solid curves represent the results of the HI 2006 DPDF Fit B to the LRG data J51 
reduced by a global factor 1.20 to correct for the contributions of proton dissociation processes. 
The dashed curves indicate the contribution of pomeron exchange in this model. 
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Figure 15: The logarithmic Q 2 derivative of the reduced diffractive cross section 
xjp a^ 3 \(3, Q 2 , xjp) at different fixed values of x P and 0. The inner error bars represent the 
statistical errors. The outer error bars indicate the statistical and systematic errors added in 
quadrature. The solid curve represents the results of the HI 2006 DPDF Fit B [O at xjp = 0.016 
reduced by a global factor 1 .20 to correct for the contributions of proton dissociation processes. 
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Figure 16: The ratio of the reduced diffractive cross section a r (f3,Q 2 } xjp) to the reduced 
inclusive cross section a r (x = f3xp,Q 2 ) obtained using the parameterisation H1PDF 2009, 
multiplied by (1 — j3)xp, shown as a function of j3 for different values of xp and Q 2 . The 
error bars indicate the statistical and systematic errors added in quadrature. The solid curves 
represent predictions obtained using the HI 2006 DPDF Fit B for the diffractive cross sections 
and HI PDF 2009 set for the inclusive cross sections. The results for the ratio derived from the 
PDF predictions are reduced by a global factor 1.20 to correct for the contribution of proton 
dissociation processes. 
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Figure 17: The ratio of the reduced diffractive cross section a r (f3,Q 2 ,xjp) to the reduced 
inclusive cross section a r (x = f3xjp,Q 2 ) obtained using the parameterisation H1PDF 2009, 
multiplied by (1 — j3)xp, shown as a function of Q 2 for different values of xp and j3. The 
error bars indicate the statistical and systematic errors added in quadrature. The solid curves 
represent predictions obtained using the HI 2006 DPDF Fit B for the diffractive cross sections 
and HI PDF 2009 set for the inclusive cross sections. The results for the ratio derived from the 
PDF predictions are reduced by a global factor 1.20 to correct for the contribution of proton 
dissociation processes. 
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Figure 18: The logarithmic Q 2 derivative of the ratio of the reduced diffractive cross section 
a? (/3, Q 2 , Xp) to the reduced inclusive cross section a r (x = (3xjp,Q 2 ) obtained using the 
parameterisation H1PDF 2009 multiplied by (1 — (3) xp, shown at different fixed values of 
xp and f3. The inner error bars represent the statistical errors. The outer error bars indicate 
the statistical and systematic errors added in quadrature. The solid curves represent predictions 
at xp = 0.016 obtained using the HI 2006 DPDF Fit B for the diffractive cross sections and 
HI PDF 2009 set for the inclusive cross sections. The results for the ratio derived from the 
PDF predictions are reduced by a global factor 1.20 to correct for the contribution of proton 
dissociation processes. 
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